Quantum trajectories: memory and continuous observation 
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Starting from a generalization of the quantum trajectory theory (based on the stochastic 
Schrodinger equation — SSE), non-Markovian models of quantum dynamics are derived. In or- 
der to describe non-Markovian effects, the approach used in this article is based on the introduction 
of random coefficients in the usual linear SSE. A major interest is that this allows to develop a 
consistent theory of quantum measurement in continuous time for these non-Markovian quantum 
trajectory models. In this context, the notions of instrument, a priori and a posteriori states are 
rigorously described. The key point is that by starting from a stochastic equation on the Hilbert 
space of the system, we are able to respect the complete positivity of the mean dynamics for the 
statistical operator and the requirements of the axioms of quantum measurement theory. The flex- 
ibility of the theory is next illustrated by a concrete physical model of a noisy oscillator where non 
Markovian effects come from random environment, coloured noises, randomness in the stimulating 
light, delay effects. The statistics of the emitted photons and the heterodyne and homodyne spectra 
are studied and we show how these quantities are sensible to the non-Markovian features of the 
system dynamics, so that, in principle, the observation and analysis of the fluorescence light could 
reveal the presence of non-Markovian effects and allow for a measure of the spectra of the noises 
affecting the system dynamics. 

PACS numbers: 42.50.Lc, 03.65.Ta, 03.65.Yz 



I. INTRODUCTION 

A first aim of the theory of open quantum systems is 
the description of the time evolution of a quantum sys- 
tem S (the open system) interacting with an environment 
£ p}. More precisely, one focuses on the reduced evolu- 
tion of S after tracing out the degrees of freedom of £. 
The resulting evolution is then usually described in terms 
of generalized master equations for the reduced density 
matrix p(t). 

A particular simple and useful way to describe an open 
system is provided by the Markovian approximation 
Essentially, this approach is based on the absence of 
memory effects in the environment. In this situation, the 
master equations are linear first-order differential equa- 
tions with a possibly time dependent generator. The gen- 
erator takes a particular form, the well known Lindblad 
form, which guarantees the complete positivity of the dy- 
namics, as required by quantum mechanics. 

Unfortunately, this approximation is no more valid 
when the memory effects of the environment cannot be 
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neglected. Several physical situations involve long mem- 
ory time effects and lead to a non-Markovian behaviour: 
strong coupling, correlation and entanglement in the ini- 
tial S-£ state 3, system at low temperature [J], struc- 
tured environment |5]. In this context, the master equa- 
tions take different forms according to the physical sit- 
uation, e.g., integro-differential equations [l|, Q, time- 
convolutionless equations [J , Lindblad rate equations 0] ■ 
The common point is that they are not in the Lindblad 
form which is characteristic of Markovian evolutions. 

In both situations (Markovian and non-Markovian) a 
useful approach to describe concrete physical evolutions 
is provided by the theory of the stochastic Schrodinger 
equation (SSE) @-[n[. A SSE is a nonlinear stochas- 
tic differential equation for a wave- function process tp(t). 
The link with the traditional master equation is given 
by the average property E[|"0(i))(0(t)|] = p(t), where E 
denotes the average over the realizations of ip(t). To find 
a SSE providing a given master equation by averaging 
is called unravelling. The idea of unravelling has been 
a real breakthrough for simulating master equations; it 
is at the root of the Monte- Carlo wave function method 
Indeed, for huge systems, the description of yj(t) 
requires much less parameters than the ones needed for 
p(t). 

However, the construction of adequate SSEs has been 
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essential also for a second aspect of the theory of open 
quantum systems: the description of the monitoring of 
S. In special situations, the SSE can be interpreted 
in terms of quantum measurements. More precisely, in 
these cases, the solution ij)(t) is called a quantum tra- 
jectory and describes the evolution of an open system 
undergoing indirect continuous measurement (continu- 
ous monitoring) [lol [l3l - [l5j |. In particular the noises 
involved in the SSE, describing jump or diffusion evolu- 
tions, can be directly connected with the outputs of mea- 
surement apparatuses. Such an interpretation is crucial 
in the understanding of real quantum optics experiments 
[Hi El EH EE Ell sucn as direct photo-detection, spectral 
photo-detection, homodyning, and heterodyning. They 
are also at the cornerstone of modern technology such as 
feedback control [g, lla [l8| . As a consequence an active 
line of research consists in finding SSEs which can be 
physically interpreted in terms of continuous monitoring 
of the system. 

In the Markovian case, this link is clearly established 
for almost all situations [13l - ll5{ . Starting from a master 
equation in Lindblad form, it is known how to construct 
an appropriate unravelling in terms of a SSE. The (non- 
linear) SSE is a stochastic equation for a random normal- 
ized vector tp(t). It is always possible to construct a linear 
SSE, driven by Poisson and Wiener noises, for a non- 
normalized vector (f>(t), such that ip(t) — <f>{t)/ \\4>{t)\\. 
Moreover, the linear and nonlinear versions of the SSE 
are related by a change of probability measure and it 
is this link which allows for a measurement interpreta- 
tion consistent with the postulates of quantum mechanics 
(lol EI| • In mathematical terms the change of measure 
is a Girsanov transformation with probability density 
||</>(i)|| 2 ; the key point which allows this transformation 
is the fact that ||0(t)|| 2 turns out to be a martingale [l4j . 
Moreover, these stochastic differential equations can be 
deduced from purely quantum evolution equations for the 
measured system coupled with a quantum environment, 
combined with a continuous monitoring of the environ- 
ment itself [ij,[20j]. 

In the non-Markovian case, to find relevant SSEs, de- 
scribing both non-Markovian quantum evolutions and 
continuous monitoring, is a tremendous challenge. Con- 
trary to the Markovian framework, no general theory has 
been developed. Essentially there exist two strategies. 

The first strategy consists in considering a physical 
model described by a non-Markovian master equation 
and in finding an appropriate pure state unravelling. 
This approach has been successfully applied in various 
situations. A major common point consists in replacing 
the memoryless white noises, used in Markovian SSE, by 
coloured gaussian noises and in introducing some delay 
effects 21]. This allows to introduce correlations in time 
which describe strong memory effects of the environment 
interaction. In this direction several models such as for 
example the so called non-Markovian quantum state dif- 
fusion have been derived [22| . Other investigations in- 
volving non-Markovian jump type SSEs have been also 



proposed (23l . |24| . While such approaches are efficient for 
simulating relevant non-Markovian evolutions, the mea- 
surement interpretation of the underlying SSEs is highly 
debated |24l - |2q ] and a complete conclusion remains still 
open. A principal problem concerns the interpretation 
of the underlying noises as outputs of continuous time 
measurements. For other models such as Lindblad rate 
equations, different type of jump unravellings have been 
proposed [27j . Jump-diffusion generalizations with mea- 
surement applications have been derived in [28[. In this 
context limitations also appear in the sense that the type 
of observables which can be measured must have partic- 
ular and restrictive forms. 

A second strategy is first to generalize directly the 
Markovian SSE by introducing memory effects. Then, 
one has to show whether this SSE provides the unrav- 
elling of some non-Markovian evolution and whether it 
has a physical measurement interpretation. To work at 
the Hilbert space level guarantees automatically the com- 
plete positivity of the evolution of the statistical opera- 
tor. In this paper we propose non-Markovian SSE models 
with physical measurement interpretation. Our strategy 
consists in adapting the Markovian approach by replacing 
white and Poisson noises by non-Markovian noises and by 
allowing for random coefficients in the equation. First 
we start with a linear SSE driven by coloured noises and 
involving random operator coefficients; the whole ran- 
domness is defined under a reference probability. Next 
we introduce the physical probability and the nonlinear 
SSE, which is a stochastic differential equation under the 
new probability. Finally it is possible to pass to the linear 
and nonlinear versions of the stochastic master equation 
(SME) and we show that they determine the dynamics 
and the continuous measurements without violating the 
axiomatic structure of quantum mechanics. The general 



mathematical structure was introduced in [29|,|30|; in 32 1 
we started to show how such a structure allows for the 
introduction of some coloured noises, while in [33| we con- 
sidered memory effects due to feedback with delay. The 
present article is devoted to the exploration and clarifica- 
tion of physical effects which can be treated within such 
a theory and to show this on a concrete physical system. 

The paper is structured as follows. Section HT1 describes 
the general theory of the stochastic Schrodinger equation. 
We present the general mathematical ingredients neces- 
sary to develop the generalization of the SSE involving 
coloured noises and random coefficients. We consider 
first a linear stochastic equation for a non normalized 
wave function <fi(t). Then, with the help of a change of 
measure, determined by the SSE itself, and a normaliza- 
tion, we derive the nonlinear SSE for the wave function 
rP(t) = (f)(t)/\\^(t)\\. In Section Hn] the linear SME and 
the nonlinear one are introduced and the measurement 
interpretation is justified by introducing positive opera- 
tor valued measures, instruments, a priori states (mean 
states) and a posteriori states (conditional states). Sec- 
tion [TV] is devoted to a concrete model, a noisy oscillator 
absorbing and emitting light, by which physical effects 
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can be discussed and the possibilities of the theory can 
be explored. Moreover, in this section we study the be- 
haviour of the outputs of the oscillator; in particular we 
study the effects of the non-Markovian terms in the dy- 
namics on the homodyne and the heterodyne spectra of 
the emitted light and on the statistics of the photons, now 
analysed by direct detection. Conclusions are presented 
in Section [V] 

II. THE STOCHASTIC SCHRODINGER 
EQUATION 

When introducing non-Markovian evolutions for a 
quantum state, the hrst problem is to guarantee the com- 
plete positivity of the evolution of the state (statistical 
operator) of the reduced system. Then, if one wants to 
introduce measurements in continuous time, the second 
problem is to have equations compatible with quantum 
measurement theory. To start from the linear version 
of the SSE allows to introduce memory by using ran- 
dom coefficients and coloured noises (no problem on com- 
plete positivity because we are working at Hilbert space 
level) and to construct the instruments related to the 
continuous monitoring (no problem with the axiomatic 
of quan tum theory because we are respecting linearity) 

[II III 111- 

Let us denote by Jf? the Hilbert space of the quantum 
system of interest, a separable complex Hilbert space, 
by jSf(Jf) the space of the bounded operators on Jf?, 
by 3T(J?) C &(Jf?) the trace class and by J^(Jf ) C 
^?(J$f) the convex set of the statistical operators. 

A. The linear SSE and the reference probability 

The starting point of the whole construction is the lin- 
ear SSE for a stochastic process <p(t) with values in Jjf: 

m 

d<f)(t) = K (t)(f)(t-)dt + L °(*)<H*-) dMj(t) 
j'=i 

d' 

+ J2(Rk(t)-t)<l>(t-)dN k (t). (1) 

k=l 

In the Markovian case the Mj are Wiener processes, 
the Nk are Poisson processes and all these processes are 
independent. Moreover, the operators Rk(t), L®(t), K(t) 
are not random. The key property, which allows for a 
measurement interpretation, is that ||0(i)|| is a mean one 
martingale and this requirement imposes a link among 
the operators R k (t), L?(t), K (t). The non-Markovian 
generalization is to take more general processes as driv- 
ing noises and to allow for random coefficients. Now we 
illustrate the precise meaning of the various quantities 
appearing in the SSE (P). 

First of all we work in a reference probability space 
(CI, J^", Q); Cl is the sample space, & the c-algebra of 



events, and Q a reference probability. The physical prob- 
ability will appear when the measuring interpretation 
will be constructed in Section Hi Bl Past and present up 
to time t are represented by the events in the sub-a- 
algebra &t C &\ the family (J?t)t>Q is a filtration of 
CT-algebras satisfying the usual hypotheses, i.e. ^ s C J?i 
for < s < t, A e & with Q(A) = implies A € J*b, and 
& t = f] T>t J? T . In (Cl, (^ t )t,Q) we have d continu- 
ous, independent, adapted, standard Wiener processes 
Bi , . . . , B d and a" adapted cadlag counting process of 
stochastic intensities ik(t) > 0, which are caglad. The 
French acronym cadlag means with trajectories continu- 
ous from the right and with limits from the left, while 
caglad means continuous from the left and with limits 
from the right. The meaning of stochastic intensity is 
given by the heuristic conditional expectation 

Eq[diV fc (t)|^i] =h(t)dt; (2) 

the stochastic intensities determine the probability law of 
the counting processes (35|. Assuming usual hypothesis, 
cadlag or caglad processes, etc., are mathematical reg- 
ularity requirements useful in a rigourous development 
of stochastic calculus; what is physically important is to 
have non anticipating (= adapted) processes. 
The m continuous processes Mj (t) are given by 

Mj(t) = f fj(s)ds + J2 f bji^dB^s), (3) 
Jo l=1 Jo 

where fj (t) , bjk (t) are complex, adapted caglad processes 
such that Vf > 0, with probability 1, J Q \fj(t) \ dt < +oo, 

Jo l^ij(^)| 2 < +oo. Some typical choices are given in 
Section [TV] 

The functions t L°(t), t i-» Ru(t), t K (t) are 
strongly caglad, bounded-operator valued adapted pro- 
cesses; to be bounded is a sufficient condition to have a 
well defined general equation [30] . For physical problems 
also the unbounded case is important and, indeed, the 
examples we shall give involve unbounded operators; the 
case involving unbounded operators, but restricted to a 
Markovian dynamics, is treated in [2(| |3l|. By allowing 
for random system operators and the general noises ([3]), 
it is possible to describe random external forces, random 
environments, coloured baths, stochastic control, adap- 
tive measurements. . . 

The SSE ([T]) is a linear stochastic differential equation 
in Ito sense. The initial condition is taken to be 

0(0) =0o with E Q [||0 O || 2 ] = 1. (4) 

Note that, by suitably choosing JPq, Q and (f>o, any 
statistical operator p € y(Jf°) can be represented as 
po — Eq[|0 o )(</>o|]- The solution <j>(t) is taken to be cadlag 
and it is unique [30j | . To write </>(£_) means to take the 
value of <p just before the possible jump at time t due to 
the counting processes. 
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By using the explicit expressions of the processes Mj, 
the linear SSE can be rewritten as 



d<j)(t) = K{t)4>(t-)dt + Li(t)</>{t-)dBi(t) 

i=l 

+ J2(Mt)-*)<t>(t-)dN k (t), (5) 

k=l 
m 

where K(t) = K (t) + ^ /j (*)£?(*) and 



and the normalized version of <fi(t) 



3=1 



Li(t)=x;^(*)6ii(*). 

3=1 



(6) 



For the physical interpretation in terms of mea- 
surements, we need ||</>(i)|| 2 to be a martingale: 

Ei 



\m)\\\&* 



- \\(/>(s)\\ 2 , for < s < t 14]. We 

shall see that this is crucial for defining physical proba- 
bilities. To this end we have to compute d ||0(i)|| 2 . Here 
and in all the formulae involving stochastic differentials, 
we have to use Ito's formula and the rules of stochastic 
calculus, which are summarized by Ito's table 

dB t (t) dBj (t) = Sij (t) dt, dN k (t)dNj (t) = S kj dt, 

dBi(t) dN k (t) = dBi(t) dt = dN k {t) dt = 0. 
Then, we get 

d 

d||</»(t)|| 2 = ^(^(i_)|(L < (t)+i i (t)t)^_))cLB i (t) 



i=l 



(#_)| (if(t)t + K{t) + Li{t) ] Li{tj) <Kt-))dt 



1=1 



J2(\\Rk(tW-)\\ 2 -H(t-)\\ 2 )dN k (t). (7) 



fc=i 



The martingale property is ensured if we have 
Eq [d||0(t)|| 2 =0. By © and Eq [dBi(t)\^ t } = 0, 
we get the restriction 

d 

K(t) = -M(t)--^L i (t?L i (t) 

i=l 
1 d ' 

+ -J>(i)(l-ifc(t) t i? fc (t)) (8) 



fc=i 



with ff(t)t = ff(t). 



B. The nonlinear SSE and the physical probability 



Let us define the quantity 



\<Kt,w)\\ ^(t.w) ifp(t,w)^0, 
z if p(t, to) = 0, 



(10) 



where z G is a non random vector with ||z|| = 1 and 
we denote by ui the generic sample point in Q, as usual. 
Moreover, we introduce the processes 

mS) ■■= 2Re<V(t_)|L i (t)V(t-)), (11) 

j fc (t) :=**(*) ||ifc(tMt-)|| 2 (12) 
By condition ©, Eq. (J7J becomes 

r d 

dp(t) =p(t_)j ^m^dS^t) 

+ I] (\\Rk(tWt-)\\ 2 - l) (diVfe(t) - ifc(t)dt) L (13) 
fc=i J 

As already said, the key property of quantum trajectory 
theory is that p(t) is a mean-one Q-martingale, which 
follows from this equation and the normalization Q of 
the initial condition [3(3, Theorem 2.4, Section 3.1]. 

a. The physical probability. Now we introduce the 
new probability measures, whose physical meaning will 
be discussed in Sect. [TTTl VA e & T , 

Pl(A) := E q \p(T)l A ] = / p(T; W )Q(dw). (14) 

Due to the martingale property of the probability density 
p(t), the probabilities P^ o are consistent, in the sense that 
P; o (f) = ^ Q (F) for F G t > s > 0. 

The new probability P^ modifies the distribution of 
the processes Bi and iVfe. A very important property is 
that a Girsanov-type theorem holds [30, Proposition 2.5, 
Remarks 2.6 and 3.5]. 

b. Girsanov transformation. Under P^ o , in the time 
interval [0, T], the processes 



Wj(t) :=B 3 (t) 



m j (s)ds, j = 



,,d, (15) 



p(t):=U(t)\\, Vi>0, 



(9) 



are independent Wiener processes, while the counting 
processes Ni, . . . , Nd> change their stochastic intensities, 
which become j\ , . . . , jd' ■ The quantities m, and j k are 
defined in Eqs. ([TT]) and ([T2]). 

Note that, if for a certain index i we have rrii(t) = 0, 
V< > 0, then Wi(t) = Bi(t); the process Bi remains a 
Wiener process also after the change of probability and 
it is independent from all the other components of W . 
For instance, from Eq. (jll[) we have mi = for all initial 
conditions when the operator \Li{t) is selfadjoint for all 
t > 0. 
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c. The nonlinear SSE. Under PT , in the time in- 
terval [0, T], the random normalized vector (fTU|) satisfies 
the stochastic differential equation 

dip(t) = K(t)ilj(t-)dt 

+ (U(t) - \ mi(t)) rl>(t-)dWi(t) 



i=l 
d' 



fe=i 



|i2fc(t)^(t_)|| 



dJVfc(i), (16) 



with ^(0) = </>o and 



1 d 



E ! ^ + 5 E(J'*(*) - i k (t)Rk(tVR k (t)). (17) 

2—1 



2 , d' 
-4- 

2 



fc=i 



To get this result one needs to compute d (l/ yjp(t) 

from Eq. () 13[) and to express this differential and dc/>(t) 
in terms of the new Wiener processes; the rigourous proof 
is given in (30j . 

At least in the Markov case, it is this equation which is 
the starting point for powerful numerical methods [H.[l2l|. 



III. THE STOCHASTIC MASTER EQUATION 

Now that we have presented the theory of the stochas- 
tic Schrodinger equation for pure states, we develop the 
analog for density matrices and we introduce the stochas- 
tic master equation (SME). 



A. The linear SME 

As in the case of the SSE, we start with a linear equa- 
tion. More precisely, from ((T|) we can derive the linear 
SME for the process a(t) := \(j>(t)){^(t)\, t > 0: 



da(t) = C(t)[a{t_)]dt + jT(Li(t)*(t-) 



i=l 



a(t_)L t (rf)dB t (t) + J2(Rk(t)v(t_)R k (rf 



fe=i 



where C(t) is the following Liouville operator: 
C(t)[r] := -i [#(*)), r] - \ £ {L^L^t), r} 

i=l 

1 d ' 

fe=i 

d d' 

+ Li{t)TLi{tf + J2 ik(t)R k (t)TRk(t? . (19) 



i=i 



fc=i 



Let us stress that this operator is random. In particular, 
this makes the solution <j{t) non Markovian since the ran- 
domness of the operator C(t) introduces a dependence on 
the past. This fact will be made explicit in the concrete 
model developed in Section UVl 

Let us note that ususal master equations (without 
the driving noises B and N), but with stochastic Li- 
ouville operators, have already been considered in the 
literature as models of non Markovian evolutions. More- 
over, these equations have been derived from unitary 
system/environment dynamics by various techniques and 
approximations, see for instance [34j . 



B. The nonlinear SME 



Note that the probability density © of P^ o with 
respect to Q can be written as pit) = Tr{cr(t)}. 
Then, we normalize a(t) by defining the state p(t) = 
a(t)/ Tr {a (t)}; when the denominator vanishes we take 
for p{t) an arbitrary state. It is then possible to show that 
pit) satisfies the nonlinear SME under the new probabil- 
ity P^ |M Remark 3.6]: 



dp(t) = C(t)\p(t-)]dt + Y{Li(t)p(t-) 
i=i 

+ p(t-)Li(tf - rm(t)p(t-j)dWi(t) 

f f R k jt)pjt.)R k jty 
+ ^VTr{i4(i)tiM«-)} P[ -> 

x(dN k (t)-j k it)dt). (20) 

Everything can be expressed in terms of density matrices 
as we can write 

rmit) = 2ReTr{L i (t)p(i_)}, 



jfc(t) = ik(t) Tr{R k itfR k it)pit-)}. 

The nonlinear SME for pit) can also be directly ob- 
tained from (1161) by remarking that 



&it-))(dN k it)-i k it)dt), (18) 



pit) = mmmi 



(21) 
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C. The a priori states and the mean evolution 

The mean state, or a priori state, is defined by 

r)(t) := Eq[a(t)] = Epr \p(t)]. (22) 

By Eqs. (JUJ and $0$) one obtains 

f,(t) = E Q [C(t) [a(t)]j = E v t q [L{t) [p{t)\] . (23) 

A major difference with the usual Markovian situation 
is that in our case this equation is not closed. In the 
Markovian case one obtains an equation of the form 
fj{t) = C(t)[r)(t)], but in our situation this is not pos- 
sible, since the operator C{t) is random and contributes 
to the mean. Formally, a closed equation can be obtained 
by using projection techniques such as Nakajima-Zwanzig 
method. This construction has been derived in [29j], but 
the final equation is essentially not tractable. 

It is then clear that the mean evolution is highly non- 
Markovian. It is important to notice that our approach 
ensures that this evolution stays completely positive. We 
then obtain a completely positive non-Markovian be- 
haviour, the memory effect being encoded into the ran- 
dom Liouville operator C(t). In particular when C(t) is 
not random, we recover the usual Markovian framework. 

D. Measurement interpretation 

In this section, we present the essential ingredients in 
order to describe the measurement interpretation of our 
theory. 

1. Observed outputs 

Let us consider a^ (i, s), i — 1, . . . , mj, rihk(t, s), h = 
1, . . . , mj which are adapted and caglad kernels and ei(t), 
Wk{t) adapted and cadlag processes. We can then define 
the following processes, which represent the outputs of 
the continuous measurement process: 

m „t 

h{t) -=J2 a tj {t, s) dMjis) + et(t), 



J h {t) := 2^ / n hk (t, s) dN k (s) + w k (t). 
k =i J (o,t] 

The idea underlying the construction of these processes 
is that the instantaneous outputs are the formal deriva- 
tives Mj(t), N k (t). The measuring apparatuses have a 
smoothing effect on the singular instantaneous outputs 
and can also provide some post-measurement process- 
ing of the outputs. These effects are represented by the 
integrals with the detector response functions aij, n^k- 



Moreover, it is possible that the detectors introduce some 
further noise, for instance of electronic origin, and this is 
taken into account by the additive noises ei, w k and by 
the fact that response functions can be random. 

Let us consider now all the events that can be observed 
up to time t, that is the events determined by the outputs 
Ig, Jh up to t. Let us call ^ the collection of such events. 
In mathematical terms is the sigma algebra gener- 
ated by h(s) and Jh(s), with s € [0, t], I = l,...,mj, 
h = l,...,mj. Because all the processes involved in 
the definition of the outputs are (^)-adapted, we get 
Sft C &ti f° r an t. Let us stress that in general we do 
not have §?< = &u because ^ contains only events which 
can be observed by the measuring apparatuses, while J^t 
can contain extra sources of noise which can affect the 
system (a noisy environment for instance). 

2. Feedback 

In this formalism we can describe also measurement 
based feedback: part of the outputs are used to control 
some features of the dynamics or of the measuring ap- 
paratus, say through a stimulating laser or through a 
local oscillator in an homo/hetero-dyne detector. When 
the feedback involves the output in the past, other mem- 
ory effects are introduced. A typical measurement based 
feedback is represented by an Hamiltonian term function- 
ally dependent on some output up to the current time; 
while in this way it becomes a random Hamiltonian, its 
contribution is perfectly compatible with the whole for- 
malism. We shall not give examples in this paper; theory 
and some applications can be found in [3CJ, Sect. 4.4] and 
& 

3. Instruments and a posteriori state 

A cornerstone of a consistent measurement interpreta- 
tion of SME relies on the introduction of the so called 
instruments. In order to develop this theory we need to 
define the propagator A(t, s) of Eq. (fT8|) . that is, the ran- 
dom linear map ct(s) i— > 5{t). An essential point is that 
this application is completely positive and satisfies the 
composition rule A(t, s) — A(t, r)oA(r, s), < s < r < t. 

Now for an event A € Sft, we define 

l t (A)[p}=E Q [l A A(t,0)[p]\- (24) 

For all A G ^t, ^t(A) is a completely positive linear map 
called instrument. In particular this gives the probabil- 
ity that an event A G Sft occurs. More precisely, if po 
represents the pre-measurement state, the probability of 
A € is given by 

TT{l t (A)[p }} = P; [A], Eq[|0o><0o|] = Po, (25) 

and we recover the previous definition of the physical 
probability. 
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Then, we can define the a a posteriori state by 
p( t ):=E P .^)|^].^, 

where 

a(t) :=E q [a(t)\& t ]. 

The state p(t) corresponds to the update of the state of 
the system conditionally to the observation of the outputs 
up to time t. 

It is important to notice that in general we can not 
derive a closed equation for p(t) such as the one for p(t). 
Essentially, it depends whether ^ t = &t or n °t- In the 
case of S?t ^ &t the randomness of the operators appear- 
ing in C{t) will prevent to close the equation; again some 
projection technique could be used to obtain a kind of 
closed equation, but intractable for practical purposes. 

As a conclusion, we can see that this approach allows 
us to describe non-Markovian evolutions which are gen- 
eralizations of the Markovian setup. As we shall see, the 
randomness of the operators Li(i), Rk(t) will be used to 
describe concrete non-Markovian effects such as coloured 
environments and incoherent stimulating light. 

The physical model is determined by the physical prob- 
ability, the nonlinear SME, the outputs, not by the SSE, 
which is not unique. Two SSEs giving solutions which 
differ only by a stochastic phase are physically equiva- 
lent; no physical consequence depends on a global phase 
in 4>{t) or Tj}(t) Sect. 2.5]. 

IV. A MODEL: A NOISY OSCILLATOR 

Let us present now a mathematically treatable, but 
sufficiently rich and physically interesting model; the aim 
is to understand what kind of physical phenomena and 
memory effects can be described by the theory we have 
presented. To be simple we take a linear system, but 
we allow for absorption, emission, coloured noises acting 
on the system and on the detection apparatuses. . . The 
general scheme is the following: 

1. The quantum system is a single oscillator; to fix 
the ideas we think of a mode in an optical cavity, 
but it could be an ion in a trap or some other sys- 
tem in the harmonic approximation. Let a and a* 
be the usual annihilation and creation operators of 
quanta in the mode; then, the free Hamiltonian of 
the oscillator is 

H (t) = H Q = v a<a, v > 0. (26) 

2. The system emits and absorbs light; the interaction 
system/electromagnetic field is treated in the usual 
Markov approximation. 

(a) Some emitted light reaches a photocounter: 
direct detection. The post-processing of the 
output is taken into account by a detector re- 
sponse function. 



(b) Some light reaches a homo/hetero-dyne detec- 
tor. The function describing the local oscilla- 
tor can be random, a way to model imperfec- 
tions. We shall show that this fact introduces 
memory in the detection process, not in the 
mean dynamics. Moreover, we can have also 
a detector response function acting as a fre- 
quency filter, see Eq. (|9ip. 

(c) We introduce a stimulating laser; the laser- 
wave/oscillator interaction is treated in the 
usual dipole and rotating wave approxima- 
tions. The laser wave can be random because 
the laser is noisy and/or because of feedback. 
This introduces memory also into the Liouvil- 
lian and in the mean dynamics, in spite of the 
fact that the interaction is without memory. 

3. We introduce various kinds of coloured environ- 
ments. According to the choices of the parameters 
these new terms can describe incoherent light, a 
squeezed reservoir, a usual (or coloured) thermal 
bath, intermediate situations. . . 

As in the general part, also in this model B\, ... , Bd 
are d independent standard Wiener processes under the 
reference probability Q; here we shall have d > 5. More- 
over, we shall introduce m = 5 diffusive channels and a 
single jump channel, d' = 1. Finally, we shall introduce a 
single diffusive output and a single counting output; ac- 
cording to the notations of Section IIII D II we shall have 
mi — mj = 1. 

A. Stimulating laser and emitted light detection 

As already said we consider the oscillator/e.m. -field 
interaction in the dipole and rotating-wave approxima- 
tions. We divide the directions of the propagating light 
into some "channels". The index 1 labels the "side" 
channels used to describe the emitted light reaching a 
photo-counter (direct detection) or a homo/hetero-dyne 
detector. Channels 2 and 3 are the "forward" channels in 
the direction of the stimulating laser; they describe also 
losses of light. 

1. Detection 

a. Direct detection. We consider only one counter, 
so that we have d' = 1. Under the reference probability 
Q, the associated counting process Ni(i) = N(t) is taken 
to be a Poisson process of intensity i\(t) = A > 0. When 
A = this channel is not open. The associated operator 
is 

R l (t) = R = ^a, /3 eC. (27) 

With respect to the general case of Section IIII D 1[ let 
us consider only deterministic, time invariant, real and 
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continuous detector response functions, so that we have 
Ji(t) = J(t), nu(t, s) = Fj(t - s), wx{t) = 0, and the 
output current is 



J(t) 



Fj(t-r) dJV(r). 



(28) 



b. Homodyne or heterodyne detection. As usual, 
homo/hetero-dyne detection is described in the Markov 
approximation by a diffusive channel driven by a Wiener 
process B\ (under the reference probability Q) [3, Sect. 
7.2]. By particularizing the quantities introduced in Sec- 
tion HE] we have L\{t) = L\(t) and M\{t) = B x (t), 
which means fi(t) = and bu(t) = Su. Then, we take 

Li{t) = -iol/i(i)a, = 1, ai G C; (29) 

hit) is the contribution of the local oscillator, which can 
be random. Randomness in the local oscillator can be 
due to imperfections, but it could be due also to the fact 
that h(t) is taken dependent on some of the observed 
outputs at previous times in order to describe adaptive 
measurements, as done in [IJ|, Sect. 7.9.2]. 

We consider again a deterministic, time invariant, real 
and continuous response function Fj ; in terms of the no- 
tation of Section fill D II we take I\(t) — I(t), a\j(t, s) — 
Fi(t — s), exit) — 0. Then, the output current of the 
homo/hetero-dyne detector is 



J(t)= / Fj(t-r)dBi(r). 
Jo 



(30) 



We assume the response function Fj to be in L 1 (M.+ ), so 
that its Fourier transform exists: 



P+OO 

Gj(M) := / e^*Fj(t)dt 
Jo 



(31) 



We shall see in Section llV El that \Gi(/j,)\ 2 has the role of 
a linear frequency filter on the output. 

c. Contribution to the linear SSE. Summarizing, 
the contributions to the right hand side of the linear SSE 
(HJ) or ([5]) of the two detection channels are 



A 



iVa) <j>(t-)dt+(0a-±) <t>(t-)dN(t) 



\ax\ t 



a* a d * - \axh{t) a4>{t-)ABx{t). 



The final linear SSE is given by Eq. ([5 



2. The forward channels 

Channels 2 and 3 represent the forward channel (the 
direction of the stimulating laser) and the lost light; we 
can include in these channels other Markovian dissipative 
contributions. There is no detector associated to these 
channels and we choose to put a diffusive component (the 



Wiener B2) in channel 2, while channel 3 is used to com- 
plete the Hamiltonian part with the contribution of the 
stimulating laser. With respect to the symbols used in 
the linear SSE (fT]) and in Section Hi AI we take 



L° 2 (t) = L° 



ioia, L%(t) = L u 3 = -ia 2 a T , (32) 



with a 2 e C, and f 2 {t) = f 3 (t) = f(t), b 2l {t) = 5 2i , 
b 3i (t) = 0, which give L 2 (t) = L° 2 , L 3 (t) = 0, 

dM 2 (t) = Jttjdt + dB 2 (t), dM 3 (t) = f(t)dt. (33) 

a. Contribution to the linear SSE. Summarizing, 
the contributions to the right hand side of the linear SSE 
([5]) of channels 2 and 3 are 

^-^y- ° ta - ii? /(*)J 0(*-)di - ma<j){t-)dB 2 {t), 

where Hf{t) contains the interaction between the stimu- 
lating external laser and the oscillator: 

H f {t) = o^W)a + a 2 f{t)al (34) 

b. Stimulating laser. The function f(t) represents 
the laser wave, eventually a laser with imperfections [lij ]. 
In the case of closed loop control, the laser wave could de- 
pend on the observed output j33| , but here we disregard 
the possibility of feedback. Then, a good model for a not 
perfectly coherent stimulating laser is the phase diffusion 
model (36}. Let v 3 > be the carrier frequency of the 
laser light (in this case Av = vq — v 3 is called detuning) 
and let e > be its bandwidth; then, 

(35) 

The quantity g contains the amplitude and the initial 
phase of the laser; in principle it could be a random vari- 
able, but for simplicity here we take it to be deterministic. 

To identify the bandwidth of the laser light /, we con- 
sider its spectrum. Being / a complex stochastic process, 
its spectrum is given by the classical definition [371 ] 



f(t) = gexp{-iiv 3 t + i^B 3 (t)} , ge 



Sf(u) := lim — Ec 



fit) At 



(36) 



By using the autocorrelation function (| A2|) of the process 
/ we get easily the Lorentzian spectrum 

SfW = ~, g l g 2 |2 — ■ (37) 



(j* - v 3 X + s 2 /4 

c. Homodyne detection. In this case the local oscil- 
lator and the stimulating light are generated by the same 
laser. A choice, which takes into account the differences 
in the optical paths, is 



h[t) 



(38) 



l/(*-A*)| ! 

we are assuming /(<) ^ 0. The phase 9 and the time 
shift At depend on the physical implementation of the 
homodyne apparatus and could be random, but, for sim- 
plicity, we take both to be deterministic, 9 € [0, 2ir) and 
At € E. 
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d. Heterodyne detection. The local oscillator h(t) 
and the stimulating wave f(t) are produced by different 
laser sources and the phase difference is not stable; the 
carrier frequencies are generally different. In this case h 
could depend on the output (another form of closed loop 
control) or could be described by a phase diffusion model 
(noise in the local oscillator). In this second case we can 
take 

h(t) =exp{i^-iz/t + iVK-B 4 (i)}, (39) 
v e R, k > 0, d e M. 

3. Summary of the contributions to the linear SME 

We have already explicitly given the various contri- 
butions to the SSE. To understand better the meaning 
of these terms it is worthwhile to write down how they 
contribute to the linear SME (fT8|) and to the random 
Liouville operator (|19l) . Let us consider the free Hamil- 
tonian of the oscillator and all the other terms we have 
introduced up to now; let us set 



7o := |ai| + |a 2 | 



A> 0, 



(40) 



and a{i) := \</>(t)) (<f>(t)\ as in Sect. 11101 Then, we have 

dcr(i) = C cm (t)[a(t-)]dt 

+ (|/?| 2 aa(t-)a* - er(t_)) (dN{t) - Xdt) 

+ (iaih(t)a(t-)(J - iwW) a ^( t -)) dB i(t) 

+ (ia 2 5-(t_)o t - ia2~aCT(i„)) dB 2 (t) + ■■■ (41) 

(the dots stay for further contributions we shall introduce 
in Section llYB)) . 



4m(«)M - -i[H +H f (t), t] 



+ 7oaTa t_^| a t flir }. (42) 

From these equations it is apparent that the elec- 
tromagnetic interaction has been treated in the usual 
Markov approximation; we see also that the parameter 
7o is the mode width. The only possible sources of mem- 
ory are / (the stimulating laser light) and h (the local 
oscillator). So, up to now the memory is due only to the 
imperfections inducing randomness in the involved lasers. 
Remember that we have not included a conceptually very 
important source of memory, the possibility of feedback. 



B. A coloured environment 

Our aim here is to introduce some sources of coloured 
noise; they could describe physically different scenarios, 
which we shall discuss at the end of the section. 



by 



Let us introduce a complex Gaussian process Y given 



Y (t) ■■= E (^w + f* x j( s ) ds ) > ( 43 ) 



3=5 



Xj(s)= / Cjis-^dBjiu). (44) 
Jo 



Here &, G C and we set 



g :=]>>A fc:=]TN 2 ; ( 4 5) 

moreover, we assume the complex functions Cj to be in- 
tegrable, i.e. 



\cj(t)\ dt < +oo. 



(46) 



Now, we add two more diffusive channels; with the 
notations of Section UH we take m = 5 and 

L°(i) = L\ = -ia, L° 5 (t) ^L Q b = -iot, (47) 



M i (t)=Y{t), M 5 {t) = Y{t), (48) 



which means f 4 (t) = f s (t) = J2j=5 x j(t), 

[h, i > 5 

This gives Li(t) = and, for i > 5, 

Li(t) = L l = blLl + hLl = -i (ha + bi^) . (49) 

The contribution of these new terms to the linear SSE 
([1]) turns out to be 



D 
~2~ 



<j)(t-)dt - i (a<dY(t) + adY(tj) 4>(t-), (50) 



where 



D := 2ka^a + k + qa 2 + qa) 2 . 



1. The spectrum of the Gaussian noise 



(51) 



The dynamics of our system involves the differential of 
the process Y or, in other terms, its generalized derivative 
Y(t). As the spectrum of / (|36|) . the spectrum of this 
classical complex process is defined by 



S Y (fi) 



lim — Ef 



: i/i *dF(t) 



(52) 
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when the limit exists. 

By construction, Y is a Gaussian process with zero 
mean; its second moments, needed in ([52]) . can be eas- 
ily computed by using the properties of the stochastic 
integrals (the Ito isometry). 

Let us introduce the Laplace transform of Cj 



Sj (z) := 



+ 00 



j(t)dt, Rez>0, (53) 



which exists due to the integrability condition (|4"6")) . 

The spectrum (|52l) is computed in Appendix IA 21 and 
it is given by 

d 

Sy(jjl) = ]T S Yj (,i), S Yj (Ai) = \bj + s 3 {-^)\ 2 . (54) 



Note that the spectrum of Y(t) is the sum of the spectra 
of the components Yj (t) without any interference among 
them. Each spectrum component contains a white-noise 
contribution (bj) and a regular one (sj), which interfere 
(they sum up inside the square modulus). Moreover, let 
us stress that by this construction it is possible to in- 
sert Gaussian noises with given spectra, not only in this 
model, but even in the general theory. 



2. Contribution to the linear SME 

As done for the electromagnetic contributions in Eq. 
(|4"Tj) , it is useful to identify the contributions to the linear 
SME due to the new noises: 



a 

da(t) = ---+Y J j£j(t)lHt-)}dt 

-i[Cj,a(t-)]dBj(t)) (55) 



i=5 



(the dots stay for the contributions already introduced), 
where 

Cj (t) [r] := -i [Hi (t), r] - ± [Cj , [Cj , t]] , (56a) 



a. Incoherent light. Consider the index j = 5 and 
assume 65 = 0. Then, C5 = and this term contributes 
only with a regular random Hamiltonian term H§ (t) . Its 
structure is very similar to Hf(t), but the two random 
processes involved are qualitatively different. The pro- 
cess / is the exponential of a Gaussian process and rep- 
resents a quasi-monochromatic wave (laser light). The 
process A5 is Gaussian and could represent, for instance, 
incoherent light with arbitrary spectrum | S5 ( — i/^) | 2 , for 
instance thermal light with a blackbody spectrum 0, 
Eqs. (1.52), (7.148)] 



M-iM)| 2 = 



Qh?^ie 



^ 2 fcjT 2 (l 



-hfi/k B T\ 



H>0. 



Another possible choice for incoherent light is an 
Ornstein-Uhlenbeck process, which means to take 



—^5 * 



75 . 
X5 = y - 1^5, 



c 5(t) = l(o,+oo)(i)ff5e 



with g 5 € C, ^5 € R, 75 > 0. In this case, from ([53|) we 
get 

S5(Z) = _ , 
Z + H5 

and the contribution to the spectrum (|54"|) is the 
Lorentzian term 



|s 5 (-iM)r 



135 1 



(^ - i*0 + 75 2 /4 



(57) 



Note that also the phase diffusion model ([35]) of the laser 
gives a Lorentzian spectrum (1371) . but, in spite of this, 
the two cases are completely different. The wave ([33)) 
is quasi-coherent, while the Ornstein-Uhlenbeck process 
represents a Gaussian incoherent wave. 

b. Squeezed reservoir. Consider now the indices j = 
6, 7 and assume cq = cj = 0; then, we get Xq — X-j — 
and the contributions of these terms are Markovian. 
Indeed, by defining 



N 2 + N 2 
70 



b 6 2 + b 7 2 
70 



(58) 



C, 



bj(J 



by 



Hj (t) := Xj (ty + Xj (t) a. (56b) 



Let us note that the contributions of the classical pro- 
cesses Xj to the dynamics (|55p are very reminiscent of the 
contribution of classical processes in the "adjoint equa- 
tion" in [l3j, Sect. 3.5]. This shows that these contri- 
butions, or at least some of them, can come out from 
the interaction of the system with a quantum reservoir 
and that could be derived by using the techniques of the 
quantum Langevin equation and the adjoint equation [l3|, 
Sects. 3.1, 3.5]. 

Now we can identify the physical meaning of various 
possible contributions. 



dC(t) 



b 6 dB 6 {t) + b 7 dB 7 (t) 



b 6 \ 2 + \b 7 \ 2 



we get 
7 



; oJa 



Cj (t) [t] = 7o« (ara) + a) ra — -{( 
j=6 ^ 



- 1! /)/ I ara — — {aa, r} 



11 



7 

= iV^([T,a}dCtt)~[a\T}dC(t)y 

By combining these contributions with the dissipative 
term in £ em (i) (|4"2"|) we get the typical dissipative effect of 
a squeezed reservoir [13|, Eq. (10.2.42)]. So, we can inter- 
pret n as the effective photon number and m as squeezing 
parameter of the reservoir. When m — 0, i.e. bj = ±i6g, 
the previous dissipative terms reduce to the contribution 
of a thermal bath with ^qu — 2 |6g| 2 . 

c. The generic case. The spectrum (|54p shows that 
the generic case is in between the two cases discussed 
above: a Markovian dissipative contribution and a ran- 
dom Hamiltonian contribution with interference between 
them. 



C. The full model 

Putting together all the contributions we have intro- 
duced in Eqs. ([2"o]l. (I5U)) and in Paragraph IIV A 2 al we 
get that the full 1SSE can be written as 

d</>(t) = ^-xoa f a+ X ~ 2 D ^J 4>(t-)dt 

- ia t 0(t_)dYi(t) - ia0(t_)dy 2 (t), 

+ (/? a -!)<£(*_) dZV(*) 5 (59) 
where A is the intensity of the Poisson process iV, 

dYi(t) = a 2 f(t)dt + dY(t), (60a) 

dY 2 {t) = cnh{F)dB 1 (t) +o^J(¥jdt 

+ o^dB 2 (t)+dY{t), (60b) 

*b = -iH> + y, IM*)I = 1; (61) 
7o, y, D, are given by Eqs. (gDJ), (031), (|ST|) . 

_Z. 27ie Liouville operator and the linear SME 

From Eqs. (|42j) and ([56]) we get the full Liouville op- 
erator 

C(t)[r] = -i[H(t), r] + 7o faro* _ I { a t aj 

1 d - - 

- 2 E [ 6 J at + 6 ^ a ' M + 6 -> a ' T ] ] ' ( 62 ) 

J'=5 

H(t) = H + f(tja + f(t)a\ 



d 

f(t) = a 2 fit)+J2 x j(t)- 

3=5 

Note that the Liouville operator is random only because 
of the presence of the process f(t), but this is enough 
to prevent from having a closed master equation for the 
a priori states, see Sect. MI CI Let us stress again the 
different physical roles of the wave / and the Gaussian 
processes Xj, as discussed in Paragraph IIV B 2 al 

By putting together Eqs. (|4ip and (|55|) we get the full 
linear SME 

d 

da(t) = £(t)[a(t-)]dt - \b i a + b i a l! , a(t)] dBi{t) 

i=5 

- i (w[h(t)aa(t) - aih(t)a(t)a^ dB 1 (t) 

— i (a<2 aa{t) — a2cr(t)a i ) d£>2(i) 
+ (|/3| 2 aa(t-)ai - ait-)) (dN(t) - Xdt). (63) 

In the measuring process there is one more source of 
memory coming in through the randomness in the local 
oscillator hit). 

2. The solution of the linear SSE 

The simplifying mathematical feature of the model we 
have constructed is that the SSE leaves invariant the co- 
herent states. Such a states are defined by ae(£) = £e(£), 
ll e (£)ll = 1- We assume the initial condition to be the 
coherent vector 

<A) = e(£o), Co e C. 
To find the solution of the SSE ([5Uf we make the ansatz 

As ||</>(i)|| 2 = \Vit)\ 2 e KcZ ^ must be a martingale, we 
ask 

• V(t) to contain only the contributions from the 
jumps of Nit) and to be such that |V(t)| is a mar- 
tingale, 

• Z(t) to contain only the contributions from the dif- 
fusive processes and to be such that e RcZ ^ is a 
martingale, 

• £ it) to be a generic stochastic process whose differ- 
ential contains all the possible terms. 

Then, we identify d(e(e)|0(t)) with (e(e)|d0(i)), where 
now d<f>it) is taken from Eq. (|59p . By equating the coef- 
ficients of the monomials with homogeneous powers of e 
we get stochastic equations for the unknown terms, which 
can be solved. The final result is 

=e-^%-iU f it)-iU Y it), (65) 
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U f {t) :=a 2 [ e-^^/^dr, 
Jo 



j=5- 



9j\ 



(t) :=e-* ot b j + / e-^-^c^uJdu, 



(66) 



(67) 



(68) 



(*) = -2i y £(s) (zsrfc(*)dBi(*) +a^dS 2 (s)) 

|(i)dyi(s) + y (ar 2 M^) 2 + c^ 2 )^( s ) 2 d s 

- (|ax| 2 + |a 2 | 2 ) fmfds, (69) 



2iRe 







l 2 l^)l 2 )d s 



re(0,t] 



AN(r) 



(70) 



Let us stress that, almost surely, the product in ([70)1 
contains a finite number of factors different from one and 
that the expression ([70| is nothing but the solution of the 
linear SDE 



dv(t) = v(t)S.^{i-m 2 \m\ 2 )dt 

+ - 1) dN(t) 



with initial condition V(0) = 1. Independently of the 
methods used to find the solution, by using stochastic 
calculus, one can check that the expression for <j)(t) de- 
fined by ([63 ]) - ([70 ]) indeed solves the SSE 

Note that the processes £(t) and Z(t) are continuous 
in time and that the contribution of the jumps is concen- 
trated in V(t). 



3. The physical probability 

Let us recall that the new probability is P^ o (dw) = 
p(t,(jj)Q(dcj) (Ti"4"l) . Having the explicit form of the solu- 
tion of the linear SSE, we can compute the probability 
density p(t) © and the normalized vector (jlOl) . which 
solves the nonlinear SSE ((15)1. From Eqs. ((9)l- ((T2)l we 
get 



p(t) = U(t)\\ =\V(t)\e 



2 RcZ(t) 



(71) 



V»(t) = exp <{ iarg^(i) + - ImZ(t) }> e(£(t)), (72) 



|V(t)| 2 =exp^ / (\~j(s))ds\ [] 

re(0,t] 



j( r ) 



AN(r) 



ReZ(t) 



2 /•* 



m. ( (s) cLB,-(s) - - rrii(s) 2 ds I , 



mi(*) = 21m \aih(t)(i(t) 



(73a) 



ma(t) = 2Im(o2^(t))- 



m=x\p\ 2 \at)f. (73b) 



Moreover, in Sect. IIVA2I we find L3 = and in Sect. 
WB\ L d (t) = and L h i > 5, ant i- self adjoint (ggj. By 
the expression (JXTJ) for mj, this gives rrii(t) = for i = 
3, . . . ,d. 

a. Girsanov transformation. As discussed in Para- 
graph III B b[ under the physical probability PT the 
process N(t) is a counting process of intensity j(t) and 
(W\,W2, -B3, . . . , Bd) is a standard Wiener process and, 
in particular, its components are independent; the first 
two components are defined by Eq. (115p . 

b. Assumption on f and h. We assume the phase 
diffusion model for the stimulating laser without feed- 
back, so that the laser wave /(<) is given by Eq. ((55)) 
and depends only on B3. Similarly we assume the local 
oscillator h(t) to be of the form (f38| or (f39)) . without 
feedback. 

c. Stochastic independence. The Girsanov transfor- 
mation of Paragraph IIV C 3 a[ the assumptions of Para- 
graph IIV C 3 bl and the results expressed by Eqs. ((65)) - 
((55)) . (f75|) have important consequences. First of all, 
under the physical probability, (WijWb) is indepen- 
dent from (B3, . . . , Bd) and £o> which is jFo-measurable, 
and, so, (Wi, W2) is independent from (£, h, /, rni,m2,j). 
Then, because (B 3 , . . . , Bd) is a Wiener process below the 
reference probability and the physical probability, the law 
of (f,h,Uf,UY,mi,m2,j) is the same below Q and P^ Q . 
Let us stress that these properties are specific of the lin- 
ear model constructed in this section, not of the general 
theory of Sections |TT] and IIII1 

d. The a priori state. From Eqs. (JUJ), ((22)), G2) and 
the properties discussed in Paragraph IIV C 3 c[ we get 
that the a priori states are given by 

V (t) = i\m)(m\] = eq [Hm))«m)\] ■ 

Note that rj(t) is a classical mixture of coherent states. 
Indeed, to keep the example simple, we introduced only 
interactions leaving invariant such a class of states. 

e. The initial condition. In the next two sections 
we shall study the detection outputs under the physi- 
cal probability PT . As we shall be interested only in the 
long time behaviour, no result will depend on the initial 
condition and from now on we take 



£0 = 0. 



(74) 
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D. Direct detection 



The mean counting intensity 



The output current of the direct detection is J(t), given 
by Eq. (|28p. where N is a counting process of stochastic 
intensity j(t) (|73b|) : 



jW = A|/3| 2 ||a^WI| 2 = A|/3| 2 |^)| 2 . 



(75) 



Let {S> t ) be the natural filtration of the processes /(•) 
and Y(-) and set @ :— Vt>o The stochastic intensity 
j(t) of the counting process N is due only to the presence 
of / and Y in £. Therefore, conditionally on fF, A is a 
time-inhomogeneous Poisson process. But, the law of 
the process j is the same under the physical probability 
and under the reference probability; this gives that the 
characteristic functional of N is 



By taking the expressions (l28l) and (|73bp of the output 
current J and of the stochastic intensity j, we get 



E pt [ J(t)] = A 



dri^(i-r)E Q |£(r)| 



As we show in Appendix IB 1 al the mean of |£(t)r has 
a limit A for large times and, by taking the response 
function Fj such that J + °° Fj(t) dt = 1, we get 



lim E r [J(t)] = lim Eq[j(*)] = A W A, (79) 



where 



^0 



= E,r 



rxp <; i / k(t)dN(t) 



exp ■ 



(e 1 ^) - l) j(t)dt 



(76) 



Similarly, one can identify all the exclusive probability 
densities. By writing 

00 pT pt n 

$$[k] = Pt(0) + V / dt n / dt n - X 
n=1 Jo Jo 

dti exp < iy^fc(tf) > pr(*n, • • -,ti), 



we have that the probability of no counts up to time T 



Pt(0) =E c 



and the probability density of a count around time t\ , . . . , 
a count around time i n and no other count in between is 



PT (t n , . . . ,ti) = Eq \j(t n ) ■ ■ ■ e" ^ 



From Eq. (|7S)) one obtains also the expressions of all 
the multi-time correlation functions of the processes N 
and J. In particular we get 



EpT 

4>0 



EpT 
*0 



k(t)dN(t) 



dt*(t)EQ [?(*)], (77) 



fci(t)dA(t) / fc 2 (s)dJV(s) 



dtfei(t)fc2(t)EQ[j(t)] 

o 

dt /" dsfci(t)fe 2 (s)E Q [?'(i)j(s)]. (78) 
o Jo 



A, 



A = A /+E A J' 

i=5 



M 2 |g| 2 (70 +f) 
70 + (A^) 2 ) ' 



(80) 



(81) 



A, = J- 



2tt 7-oo l^o+i^T 
Sj is the Laplace transform ([53 



a. Incoherent light. In the case of incoherent light 
with Lorcntzian spectrum (|57[) . described in Paragraph 
IIV B 2 al the contribution to A is 



Lggj (70 + 75) 

7075 |xo + >?5| 2 



(83) 



note that it has same form of the contribution of /. 

b. Squeezed reservoir. In the case of the Markovian 
squeezed reservoir of Paragraph flV B 2 bl we get s 6 = S7 = 
and 



A 6 + A 7 = n, 



(84) 



where the parameter n is defined in Eq. (|58l) . Note that 
the squeezing parameter m does not contribute to A. 



2. The Mandel Q-parameter 

To compute the full statistics of the counts, say all 
the exclusive probability densities, is not easy. A simple, 
significant parameter related to a counting statistics is 
the Mandel Q-parameter defined, for to > and t > 0, 
by 



Qto(t) 



Var pt+ * [N(t + t )-N(t )] 
E pt+t0 [N(t + t ) - N(t )] 



(85) 
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As for a Poisson process this parameter is zero, in quan- 
tum optics it is usual to say that in the case of a posi- 
tive Q parameter one has super- Poissonian light and sub- 
Poissonian light in the other case. Sub-Poissonian light 
is considered an indication of non-classical effects. By 
Eqs. dill), (Hi) we get 



Qt (t) 



Var Q Jt* +t ° i(-s)ds 



E; 



> 0. 



The emitted light is always super-Poissonian; indeed, it 
is well established that non-classical light can not be ob- 
tained by a linear system. 
For large times we can define 



Q(t):= lim Q t0 (t). 

to— ¥+00 



(87) 



This parameter is studied in Appendix IB 1 b| its general 
expression is given by Eq. (IB2|) . We have that Q(t) is non 
decreasing, Q(0) = and that the limit lim t _^ +00 Q(t) 
exists. The expression (|B2[) for Q(t) shows that the con- 
tributions of / and Y interfere and the same holds for the 
various components of the process Y: the situation for 
the Mandel Q-parameter is much more complex than for 
the mean intensity A, to which each channel contributes 
independently from the others, as seen in Section HV D II 



E. Heterodyne/homodyne detection 

In terms of the new Wiener process, the hetero/homo- 
dyne current I(t) ([50)1 can be written as 



I{t)= Fi(t-r)[mi(r)dr + dWi(r)], (88a) 
Jo 



mi(t) = 2 Re (-ialh(t){ip(t)\a>ip(t)}) 

= 2Im(W/i(t)£(i)) . (88b) 



As discussed in Paragraph IIV C 3c\ in our linear model 
the processes Wi and mi turn out to be stochastically in- 
dependent under the physical probability, so that, mod- 
ulo the detector response function Fi, we can interpret 
the output (|88ap as signal, mi, plus independent white 
noise, W\. 

By using Eqs. ([65 ) -pg )) . ([5"8]) and the fact that U Y 
has zero mean and it is independent of h, we get the 
expression of the mean output current 

E r t [/(*)] = -2 Rear a 2 / ds f drFAt-s) 
f° Jo Jo 

xe-^ s -^E Q [h(s)f(r)]. (89) 



Let us study now the spectrum of the hetero/homo- 
dyne output, which can be obtained experimentally by a 



spectrum analyzer. The current I(t) (188a[) is a classical 
stochastic process and again its spectrum is given by the 
classical definition 37] 



Siifi) = Km Tf^vj, 



jfit 



I(t)dt 



(gg) In Appendix IB 21 we show that we have 



S i (v) = \G i (li)\ 2 [1 + SM}, 



(90) 



(91) 



where Gi (//) is the Fourier transform (|3Tj) of the detector 
response function -Fj(t), 1 is the constant contribution of 
the white noise W% and 



T-++00 1 



mi 



{ty^dt 



(92) 



From the expression (|91j) we see that the response func- 
tion Fi acts as a frequency filter |G/(/z)| ; for instance, 
we can take any band-pass filter. Moreover, we have al- 
ways Si(jjl) > \Gi(fx)\ . This is due to the absence of 
correlations between W\ and mi and it is interpreted 
as absence of squeezing in the quadratures of the emit- 
ted light, in agreement with the fact that linear systems 
cannot generate non-classical light. On the contrary, a 
nonlinear system, such as a two-level atom, can generate 
squeezed fluorescence light and this can be treated by the 
formalism of SSE and continuous measurements [Til . [33| . 

As discussed in Appendix lB 21 the spectrum S m (ii) can 
be computed and we obtain 

S m (n) = Sii(n) + 5i 2 (m) + S 2 {v4 + fJ.) + ^2(^4- m)> (93) 



where the first two contributions are given by Eqs. (|B3I) . 
(|B4|) and 



s 2 (aO = 



+ °° , 2 74 |ai| 2 E J i 5 l^ + ^(-^)| 2 
-00 7T (7 4 2 + 4 (fj, - x) 2 J |x + \x\ 2 



(94) 



here 1/4 = v and 74 = k for heterodyning or = ^3 and 
74 = s for homodyning. 

Let us stress that S\\(p) + S\<x{}i) is the contribution to 
the fluorescence light of the laser wave f(t), while 62(^4 + 
A*) +S2 {va — A 4 ) i s the contribution of the environment and 
of the incoherent light. Moreover, by comparing Eqs. 
and (19411 . we obtain 



1 

47T 



+00 



S 2 {y)dn 



ad 



d 

3=5 



(95) 



1. Heterodyne spectrum 

In the case of heterodyning the processes / (|3"5j) and h 
are independent and this simplifies the computations 
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of the fourth order moments involved in the expressions 
(|B3|) . ([B4]) . In Appendix [ETal we show that 



S i2 (ji) =0, 



Sn(fi) = Si(v + n) + S-i(v- /i), 



(96) 



(97) 



Si(v) = \ai\ 2 \a 2 \ 2 \g\ 2 



1 



e 



(7Q+e) 2 



(Ai/)' 



(70+k) 2 



+ (l/ - i/ ) Z 



(^ + (A,) 2 ) (H^ + (---o) 2 ) 

_j_ 7o + K + f 

70 M + fa - ^) 2 



(98) 



Recall that = — ^3 is the detuning. 

By explicit computations, we find also that the mean 
current (|89|) does not contribute to the spectrum: 



lim — 

T-S-+OC T 



e^'EpT [I(t)]dt 



00 



0. 



a. Intensity. Note that from Eqs. (|8Tj) and (|98|) we 



get 



1 /" +0 ° 



(99) 



Equations ([95]). (|99|) connect the intensities in direct de- 
tection and heterodyning, a general relation alread y en - 
countered in other systems as a two- level atom [iH l3Slj - 
b. Perfect local oscillator. To summarize, the het- 
erodyne spectrum is given by Eqs. (|9"T|). (]9"3"]) , ((Ml) . (|M|. 
(]9"T]) , ([95]) with 1^4 = v and 74 = k. These formulae are 
a little bit involved, as they contain contributions from 
various sources of noise: 70 is the natural with of the 
cavity mode, e is the width of the stimulating laser, k is 
the width of the local oscillator. The situation becomes 
more transparent in the limit of a perfect local oscillator: 
k I 0. By recalling that 



— r, T2 — 777T — > 8{x + v), in 



this limit the heterodyne spectrum reduces to 



S m (p) = ^ [Si( u + +Si(v- fi)] , 



(100a) 



S 1 {v±n) = ]ai] ' Q 2 al 5/(1/ ± M ), (100b) 

(vip-m) +7o/4 



the spectra of the stimulating laser wave / and of the in- 
coherent noise Y(t) are given by Eqs. (|37|) and (f54| . Let 
us recall that Vq is the resonance frequency of our quan- 
tum system and v is the frequency of the local oscillator 
in the detection apparatus, which can be adjusted by the 
experimenter. We can say that from the heterodyne out- 
put we can see part of the input spectra Sf, Sy through 
a window of width 70 centered on vq. In principle, in the 
case of a bad cavity with a big width 70, we can read 
the Markovian and non Markovian character of the var- 
ious contributions to the dynamics from the heterodyne 
spectrum. 



2. Homodyne spectrum 

The homodyne spectrum is given again by the general 
formulae ([90]) , (]9"T1) , (192]) . Now, in the homodyne scheme, 
the local oscillator and the stimulating laser wave come 
out from the same source, in order to maintain phase co- 
herence during the detection process and to get a phase 
sensitive measurement procedure. However, this does not 
give coherence between the local oscillator and the vari- 
ous sources of noise described by the process Y{t) (ther- 
mal reservoirs, incoherent light. . . ) and, indeed, the re- 
lated contribution S^O- 4 ) to the homodyne spectrum is 
the same as in heterodyning and it is given by Eq. (]94[) 
with the substitutions 1/4 — > v 3 and 74 — \ e. 

The situation is different in the case of the compo- 
nents Su(fJ>) (]B3|) and Si2((i) (]B4[) . which contain both 
the stimulating laser wave / (]3"5"]l and the local oscillator 
(|38j) . To simplify the computations we consider only two 
limiting cases. The first case is when the optical paths of 
laser wave and local oscillator are much larger than the 
coherence length of the source; this is done by taking a 
large time delay At and this washes out any phase sen- 
sitivity. The second case is when the two optical paths 
are perfectly balanced (At = 0) and the phase sensitiv- 
ity is maximal. The general case is in between these two 
extremes. 

a. The case At — > ±00. By analysing the expres- 
sions (IB3|) and (JB4|, one can see that in this limit / 
and h become stochastically independent and Su and 
S12 become the same as in the heterodyne case with 
the substitutions v — > v% and n — > e, i.e. Si 2 (fj) = 0, 
Sxi{n) = Si(is 3 + fi) + Si(v 3 - fi) and 



+ 



g|«l| 2 l«2| 2 |ff| 2 

e 2 + fi 2 



(70 +e) 



(7Q+e) 2 



{Aut^Y 



(Auy 
1_ 70 + 2e 
70 e 2 + [i 2 



e 2 |ai| 2 \a 2 \ 2 \g\ 2 



S 2 {v ± fx) 



\<*i 



(is±fi-is ) 2 + lo 2 /4 



S Y {v±n)\ (100c) (Ilo+ili + (A^) 2 ) ((24£li + (A^ T /i) 2 ) 



(101) 
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b. The case At = 0. As shown in Appendix IB 2 bl 
in this case we have 



random Liouville operator appears and this prevents to 
have a simple closed equation for the mean dynamics. 



(Ai/) 



- Re 

7o 



70+e 



i(Au-n) 3d£-i(Ai/ + M) 



x 1- 



7oe 2i < 



7o + 2e - 2iAi/ 



, (102) 



By constructing positive operator valued measures and 
instruments, we have also shown that our proposal is 
compatible with an interpretation in terms of continu- 
ous monitoring of the system; the axiomatic structure of 
quantum mechanics is respected. The key point in the 
construction is that the starting point is the linear SSE 
and that its structure is such that the square norm of the 
solution is a martingale. 



arg 



70 +£ 
2 



iA^ 



(103) 



Let us stress that At = means that the two optical 
paths are perfectly balanced and homodyning allows for 
perfect interference, which generates a delta-spike at zero 
frequency. By varying |At| from zero to infinity, one goes 
from (fTU2f to (fTUTj) and Si 2 (p) = 0. 

From Eq. (|102[) we see that, by adjusting the phase £, 
we can change the relative intensity of the delta-spyke 
and the regular part. To simplify, let us consider only 
two extreme cases. We take always a vanishing detuning, 
which implies, in particular, £ = arg (ot\ S2) + 9. First 
we define 



16 |ai| 2 |a 2 | 2 \g\ 2 



To 



Then, for Av — and £ = ±7r/2, we have 



Snip) + Snip) 



7o + 2e 



l{p). 



In the other extreme case, Av — and £ = or %, we 
have 



Su(p) + Suip) = 



e 2 Kn) 



(70 + e) (70 + 2e) 2 



o7oi(0)<J(/x). 



Let us stress that this phase sensitivity is the character- 
izing feature of the homodyne detection. 



V. CONCLUSIONS 

In this article we have presented a SSE in which the 
the involved operators are allowed to be random and to 
depend on the past; moreover, the driving noises can 
be coloured diffusive processes and general counting pro- 
cesses, not only white noise and Poisson processes. This 
modification introduces memory without violating the 
complete positivity of the dynamics for the reduced state 
(the a priori state or mean state) . In this way we have an 
unravelling of a CP dynamics with memory. The main 
difference with respect the Markovian case is that now a 



On the physical ground, this theory allows to introduce 
in a consistent way various memory effects. The most im- 
portant one is the possibility of introducing measurement 
based feedback, with delay, a fact which opens the way 
to a general treatment of quantum closed loop control. 
Memory and feedback can be introduced also in the de- 
tection part (say in the local oscillator of a hetero/homo- 
dyne detector) and this allows for a consistent treatment 
of adaptive measurements. 

In order to understand the possibilities of the the- 
ory, we have introduced the simplest quantum system, 
an harmonic oscillator with a dynamics leaving invari- 
ant the coherent states. The possibility of using random 
terms in the Liouville operator allows to introduce vari- 
ous memory effects, such as coloured environments, pro- 
ducing non-white thermal-like effects with any spectral 
density. Moreover, the stimulating light can be a laser 
with imperfections (we have taken the phase diffusion 
model of a laser wave) and/or incoherent light (thermal 
light, for instance). 

As possible continuous monitoring, we have first stud- 
ied the direct detection and we have shown how the de- 
tected intensity and the Mandel Q-parameter depend on 
the characteristics of the input (light and thermal ef- 
fects). 

Then, we have studied the spectrum of hetero/homo- 
dyne current. Here, a new point is that also the imper- 
fections in the local oscillator can be consistently intro- 
duced in the theory. The two spectra are affected by 
all the noises and characteristics of the dynamics and 
have a somewhat complicated expressions; however, two 
features appear. The heterodyne spectrum reproduces 
a part of the spectra of the thermal noises and of the 
stimulating light and, so, in principle, it allows to see di- 
rectly the sources of non-Markovian effects. The homo- 
dyne spectrum turns out to be phase sensible, as known; 
by controlling an initial phase it is possible to change 
the field quadrature which is detected. The possibility 
of introducing a coherence length in the stimulating light 
and in the local oscillator allows to study the influence of 
imperfections on homodyning; one sees that, as soon the 
apparatus is not perfectly balanced, the homodyne spec- 
trum becomes the same as the heterodyne spectrum. 
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Appendix A: Some autocorrelation functions 
1. The stimulating laser 



Let us consider the process (|35|) . modelling the stim- 
ulating laser light. We can say that / is a log-normal 
process, i.e. the exponential of a (complex) Gaussian 
process. In stochastic calculus it is well known that 
exp {iy/e B^{t) +et/2} is a mean one martingale, which 
gives the mean function 

% [/(«)] - se-^ 4 "^ 2 



and the autocorrelation functions 



mm =\g\ 



— \n\ 2 „if3(s-r)-e|s-r|/2 



(Al) 



(A2) 



e q [mm] = 9 



_ 2 -(iws+e/2)(«+r)-e(Mr) 



(A3) 

where a A 6 is the minimum between a and b. In the 
proof of Eqs. (|A2|) and (|A3|) one has to use (IA1|) and the 
independence of the increments of the Wiener process. 

Let us stress that from (|A2|) we get immediately the 
spectrum (|3"?| . while from (|A1[) and (IA3|) we get 



lim — 

T^+oo T 



lim — Er 

T->+oo T 1 



[/(*)] d* 



e^f(t)dt 



0, 



We shall need also the autocorrelation function of the 
related process Uf(t) for large times. From (pu) . (|A2I) . 
131) we obtain 

lim Eq [Uf{t + t )Uf(s + t )} = 0, (A4) 

t — >+oo 



Then, by using the Ito isometry we get 

m 2" 



e^dY^) 



f T dt b j+ fdsc^e^ 

1 j=5 Jo J ° 



By the limit T — > +oo we get immediately the Laplace 
transform (|53|) inside the square modulus and (|54j) is 
proved. 

The expression (|43|) of the process Y contains the func- 
tion Cj(t) and its spectrum involves the Laplace trans- 
form Sj(z) (|53p of Cj. In the following, we shall need also 
the Fourier inversion formula: 



Cj(t) 



1 

2^ 



"^•(-i^d/i, t>0. (A6) 



Similarly to ([55]). (|A6|) . for the function £,(t) JBSJ we 
have 



e-\(()d< = 



*t„ ^ J* _ 6 i + *j( a; ) 



Rez>0, (A7) 



For the related process J7y(t), from Eqs. (JBTJ) , (|AT| . 
()A8|) we obtain, with t > s, 



lim 1 

to— S- + 00 



U Y (t + t )U Y (s + t ) 

d 



1 / >+ °° _W+_ S ) \ "» |&j + 

2tt 



j=5 



+ 1X1 



■ dx, (A9) 



lim E c 



Uf(t + t Q )U f (s + t ) 
a 2 \ \g\ I e 



_,— ie 3 (t—s) 



e -*3(t-s) _ e -«o(t-«) 



>*0 - 



f > s, (A5) 



where X3 = 



ii/ 3 . 



lim Eo[C/y(t + to)t/y(s + to)] = tt" / dx 
tn^+oo 2?r /_ 00 

-te(*- s ) + Sj(-ix)} [bj +Sj(ix)} 



i=5 



xo 2 + x 2 



(A10) 



3. The local oscillator in the heterodyne scheme 



2. Autocorrelation functions of the process Y 

Let us consider the Gaussian process Y defined by Eqs. 
(|33j): we can write 



/ e iAlt dY(t) = V / dB j (i)e i,li 
Jo j=5 Jo 



j /,,■ + [ T *d S e^ Cj ( S ) 



For the moments of the local oscillator in the hetero- 
dyne measurement scheme Q39p we have formulae analo- 
gous to those of Section |A~D 



E Q [h(t)} 



Eq h(t)h(s) = e 



_ Jv(t-s)- K ,\t-s\/2 



Eq h{t)h{s) = e~ 2l1> e 



— 2iif (ii/-K/2)(t+s)-«(tAs) 



(All) 



(A12) 



(A13) 
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4. The local oscillator in the homodyne scheme 

In the homodyne measurement scheme the local oscil- 
lator h is proportional to the shifted stimulating laser / 
and it is given by Eq. i.e. 



h(t)h(s) 



_,ii/ 3 (t-s)-e|t-s|/2 



(A16) 



h(t) = ¥ e 



9 e -ii^(t-At)+iVsB 3 (t-At) 

Is 

As before we get 

Eq [h(t)] = 



(A14) 



Er 



h{t) h{s) 



9 e \ (11/3-f )(t+s-2Ai)-e:(tAs-Ai) 

(A17) 



9 e -ii/ 3 (t-At)-e|t-At|/2 
i i e , 



In the homodyne case we need also the two fourth order 
(A15) moments 



J 



E r 



/(ri)/(r 2 )/i(t)Ma)] = l»| 2 ex P { il/ 3 (* - n - s + r a ) - | <5(n,t - At, s - At,r 2 )} , 
f(n)f(r 2 )h(t) ft(s)] = | ff | 2 exp {ii/ 3 (t - n + s - r 2 - 2At) - | *(n, i - At, r 2 , s - At)} , 



(A18) 
(A19) 



where 



e X p{-|S(tl,* 2 ,* 3 ,* 4 )} =Eq [exp{i^( J B 3 (t 1 )-S3(t2) + S3(t3)-S 3 (t4))}] • 

By some long, but straightforward computations we get 



<5(tl,t 2 ,t 3 ,t 4 ) 



1*4 


-* 3 |" 


f 1*2 -tl|, 


for ti V t 2 < *3 A t 4 or t 3 V t 4 < ti A t 2 , 


1*4 


-*l|- 


H*3-*2|, 


for ti V t 4 < t 2 A t 3 or t 2 V t 3 < ti A t 4 , 


1*4 


-*l|- 


f3|t 3 -t 2 |, 


for t 4 < t 2 < t 3 < ti or ti < t 3 < t 2 < t 4 , 


1*4 


-* 3 |- 


F3|t 2 -ti|, 


for t 4 < t 2 < ti < t 3 or t 3 < ti < t 2 < t 4 , 


1*2 


-*3|- 


f- 3 |t 4 — *1 1 , 


for t 2 < t 4 < ti < t 3 or t 3 < ti < t 4 < t 2 , 


1*2 


-*l|- 


f 3 |*3 - t4 , 


for t 2 < t 4 < t3 < ti or ti < t3 < t 4 < t 2 . 



Appendix B: Detection 
1. Direct detection 

a. The mean counting intensity 

The expression of £(t) is given by (1651) ; moreover, the 
processes Uf (|66|) and C/y- ([67]) are independent and Uy 
is Gaussian with mean zero. So, we have 

EQlleWlVEQll^WI^+EQllC/yWI 2 ]. 



b. The Mandel parameter 



By using Eqs. {JHSJ) , (EZJ), CMI and by recalling 

that the process Uy is Gaussian with zero mean and 
Uy(s) and Uf(r) are independent, we get 



Q{t) = 



2A|/3| 2 7 o 
At 



From (|A5|) . (|A9l) we get 



rt+to 

x lim / ds / dr 

to— >+oo 



to 



to 



|a2| 2 |ff| 2 (7o + e) 



|E Q [C/y( S )C/ F (r)] 



Covr 



CV(a)J7y(r) 



|?7/( S )| 2 ,|C//(r)r 



To 



(70+gr 



— y 



i=5 



(A^) : 



I i ■ j 2 
-oo | ^o + ^l 



■da:. (Bl) 



From Eq. (|B1[) we obtain the expression of A (|80l) . 



2RcE Q Uf(s)U f (r) E Q J7y(s)J7y(r) 



We see that Q(0) = 0; moreover, for large times the 
integrand depends only on the difference s — r and one 
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can check that Q(t) > 0. So, Q(t) is non decreasing and integration variables, 



x lim 

to— v+oo 



U Y (t + s)U Y (t ) 



\E Q [U Y (t Q + s)U Y (t ) 



2ReE r 



U f {t + s)Uf(t ) E Q U Y (t + s)U Y (to) 



\U f (t + s )\ 2 ,\U f (t )\< 



(B2) 



By us ing ([Ml) , (|X9)l . (|A10|) and the analogs of (|A18|) and 
(|A19|) . one sees that lim t _j. +00 Q(t) exists and it would 
be possible to compute Q(t), but its general expression 
is very involved and not particularly instructive. 



2. Homo/hetero-dyne spectra 



By (|88ap and the independence of W± and mi under 
the physical probability PT , we get, by some changes of 



— EmT 

T r fo 



e ltlt I(t)dt 



o 



dte ilMt F{t) 







+ -E 
T 



dte^m^t) / dse i ^F(s) 
Jo 



Then, in the limit T +oo, we obtain Eqs. flU), (l9~2"1) . 



To compute the spectrum S m (fi) (|92p we have to take 
the expression (|73a|) of mi and to use Eqs. ([65|) - (|68|) . By 
recalling that U Y is Gaussian, with zero mean and inde- 
pendent from h and /, by some computations we arrive to 
the expression S m (m) = Sn(^) + Si 2 (^) + S 2 i(^) + S 22 (fi) 
with 



5u(/i)= T lim ilf^llJM. Re / dt [ d* ••• >r "^ / «lr, / dr. <■"■" J + "" r ' 



o Jo 
t 



o Jo 



f(n)f(r 2 )h(t)h(s) + ( I dn / dr 2 + / dn / dr 2 ) e^+^Eq f (n) f (r 2 ) h(t) h(s) }, (B3) 

s JO JO JO 



4alW /" T /■* / /"* f s 

S12 (m) = nm R- e m / d M ds cos ut (t — s)) I / dri / dr 2 



+ 2 / dri ; ur 2 
'o Jo 



2 / dn / dr 2 c- >f0 ( t+s - ri - I ' 2 )E r 



f(n)f(r 2 )h(t) h(s) , (B4) 



S 2 i(n)= lim 4 ^ VRe / dt / ds / dr cos (/i (t — s)) <7,(i — r)gj(s — r) Eq h(t)h(s) , (B5) 



T^+oo T — .;„ .,„ 



T pt r s 



m 2 — '■' 



3=5 



JO JO 



S 2 2(»= lim Re ^-y^ / dt ds dr cos (fi(t- s))gj(t-r)gj(s-r)E Q h{t)h{s) . (B6) 



Let us consider now both heterodyning and homodyning; the moments (|A12j) and (|A16[) can be written in a unified 
way as 

Eq [ft(i)/l(s)] = e -* 4(t - s \ t > s, 
with 2*4 = t| — iv for heterodyning, x 4 = | — iv 3 for homodyning. Then, from (IB5[) we obtain 

SM =4| ai | 2 VRe / dt / d S e-^(*- s )cos(/i(t- S )) 5j (t)^). (B7) 
J=5 Jo Jo 
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By inserting (IA8|) into (|B7|) and, then, by using (IA7|) for gj(s), we get §2i(p) = 5*2(^4 + /x) + ^2(^4 — /x) with 52 (/i) 
given by (|94|) . 

By similar computations, from (|B6[) we find that 622 (/^) vanishes for both the choices (|38|) and (p9|) for the local 
oscillator this ends the proof of (|93p . (IM1) . 



a. Heterodyne detection 
By using (|J2|) and (jA"l"2l) inside (|B3]), we get (JW]) with 



- i T * jf - = < "- i)< - > { (/' - r *> + 1: * 1: 



dr* 



x e -xo(i-ri)-"o(s-''2)+i>'3(r2-ri)-5(ri-r2) | / / fa. g-xo(t-r2)->«o(s— ri)+ii>3(ri-r2)- f (ri -r 2 ) 

./o Jo 



= Re 



2M 2 |a 2 | 2 M 2 



1 



w 3 



1 



2 



1 



+ i (1/3 - z^) 70 



1 



( xq + § - i^3 J 7o 



By explicit computations, from (|B8)) we get ([Ml)- Similarly, by (fB4t . (|A"3|) . (|A13|) we obtain Si 2 {p) = 0. 



(B8) 



b. Homodymng, At — 
By inserting the expression (|A18|) into (|B3I) and (1A19[) into (|B4I) . we obtain 
o / \ K| 2 K| 2 |g| 2 f. e(7o + e)r 1 



To 



(7Q+g) 2 



+ O + Az/) 



2 + (70+e) 2 



(/x - A^ 



(B9) 



S 12 (,) = Re j g|2e2iVW J 4^)- ^ 



(x + w 3 + §, 



xo + 1^3 + ' 



1 



1 



>*o + i (i>3 + M) + I x o + i (^3 - M) + 



(BIO) 



In doing the computations we have used the well known 2( ^ s '"^J/ 2)2 = 1 T ™f T -> <$(zx). By adding the expres- 
approximations of the Dirac delta s -^f- -> and s i ns ([B9} and (iMil) we get (1T021) . 
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